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Abstract 

Quantum mechanics in a one-parameter family of volcano potentials is investigated. After a 
discussion on their construction and classical mechanics, we obtain exact, normalisable bound 
states for specific values of the energy. The nature of the wave functions and probability densities, 
as well as some curious features of the solutions are highlighted. 
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Introduction : Volcano potentials are not quite common in traditional quantum mechanics. A 
generic potential of this kind has a depression (well) at the origin with its value approaching 
negative infinity asymptotically (on both sides of the origin). The dip at the origin with 
finite-height barriers on either side is the reason behind the name volcano jl|, Q]. One 
might have another type of volcano potential for which the asymptotics is different-it goes 
to zero, or, a constant value, asymptotically. The former type is the inverted version of 
the well-known double well whereas the latter is the so-called double barrier In this 
article, we focus on the former type of potentials, i.e. the ones which go to negative infinity 
asymptotically. It must be noted that a crucial difference between the two types is the fact 
that the former, by definition, has non-Hermitian boundary conditions j^. 

Recently, volcano potentials have arisen in various contexts. Both the types mentioned in 
the previous paragraph appear in high energy physics in the context of localization of spin 
zero, half and spin two fields in the currently popular braneworld models For example, 
in the five dimensional Randall-Sundrum model, the equation for the Kaluza-Klein modes 
of the graviton reduces to a Schrodinger-like equation in a volcano potential of either type 
(using a coordinate transformation one can go from one type to the other, though this 
may not be possible always). On the other hand, double barrier structures (volcano box 
potentials) arise in studies related to artificial quantum heterostructures (quantum wires, 
wells and dots) and are reasonably well-known [3]. It is therefore relevant to look at model 
volcano potentials and try to understand standard quantum mechanics in their presence. 
Several authors have made such attempts in the recent past ll, 121]. Our aim here is to provide 
further illustrations and insight for a specific one-parameter family of volcano potentials of 
the first type, primarily through a class of exact solutions. 

What do we expect quantum-mechanically, if, say, a quantum particle feels such a po- 
tential ? The natural answer is that resonances appear in the spectrum which correspond to 
complex eigenvalues. The resonant states have finite masses and widths. They may tunnel 
out and disappear - but using known formalism we can estimate how long they might exist 
and provide the illusion of a possible bound state. Extensive work on understanding the 
resonances of the volcano have been reported in On the other hand, it is also not 

impossible to have bound states with a real spectrum. Examples such as the well known 
ax'^ — hx'^ potential have been studied quite extensively in the context of the recently formu- 
lated PT symmetric quantum mechanics where the non-Hermiticity of the Schrodinger 
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operator gives rise to distinct features not quite known in standard quantum mechanics. 
The volcano potentials: We begin by proposing a somewhat general construction of a volcano 
potential using arbitrary but well-behaved functions f(x) and h(x). One can obtain various 
profiles for the potential choosing different forms of these two functions. In order to generate 
a volcano potential we must have specific properties for f(x) and h(x). We delineate these 
below for symmetric potentials. 

(i) f(x) is an even function, is finite (may be zero) at the origin {x = 0) and increases 
monotonically (to positive infinity) ±00. 

(ii) h(x) is also an even function, is finite (may be zero) at the origin {x = 0) but decreases 
monotonically (to negative infinity) as x — ±00 at a rate faster than the increase of f{x). 

As an example, let us consider f{x) = ax"^ and h{x) = —bx"^, where a > and b > are 
two constants. The resultant V{x) = f{x) + h{x) is the inverted double well and represents 
a volcano potential. It is easy to construct many other examples. 

A couple of points about classical mechanics in this system. 

(i) A particle with energy E < Vmax and with initial position within —xi < x < xi 
(xi, —Xi and X2, —X2 are the four points of intersection of the E = constant line with the 
potential, Xi < X2 and E = V{±Xi,±X2)) will oscillate about x = 0. It will reside in 
the volcano and cannot come out. Within the barriers, i.e. for initial conditions between 
Xi < X < X2 OT —Xi > X > —X2 there are no classical solutions. However if the initial 
position is such that x > X2 or x < —X2 then the particle will disappear to positive or 
negative infinity. A nice discussion on the phase fiow, fixed points and invariant sets of a 
first order autonomous system in a generic volcano potential of the type — x'^ is available 
in Q. 

(ii) If E = Vmax then the particle resides at an unstable point and depending on the 
initial velocity can either slide off to ±00 or down into the volcano. The usual 'Euclidean 
time' instanton solutions for the double well will now be valid for 'real' time. 

Another way to construct a family of potentials which includes the volcano, is by using 
a single function and its derivatives. Let us write down an expression for such potentials 
using the following ansatz: 



where, ai and 02 are two arbitrary constants. We assume ^'(x) to be such that e^^*^^^ and 




(1) 
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g"{x) are positive for all x. For diverse choices of the function, g(x) one can obtain different 
potentials, each of which will have a definite shape (single barrier, volcano, single well and 
double well). g{x) therefore, is some kind of generating function for the potentials. For 

example, let us consider g(x) = In (sech x). In this case, the potential has single barrier and 
is proportional to sech^x. Similarly for g(x) = ln(cosh'' x), the general form of the potential 
becomes 

V{x) — — (oi cosh^'^ x + a2sech^a;) (2) 

where, u is an arbitrary constant. Now, for distinct values of the constants the potential 
takes various forms. Let us discuss some of the possibilities in more detail. 
Case I: ( and 02 both positive ) 
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FIG. 1: The potential, 2V is plotted as a function of x for different values of the parameter for 
ai = 0.009 and 02 = + 1). The depth at x = increases with increasing v. 

In this case one can obtain a double-barrier volcano potential with V(x) negative for 
all values of x. It is possible to have bound states at the location of the well. Assuming 
a2 — v^u + l) we find a minimum at x = if 02 > vai. The profile depends on the parameter 
V [Fig. (1)]. For the particle energies lesser than the depth of the potential well, there will 
be bound states. The potential decreases again beyond the turning points. Thus, there 
is a finite probability of tunneling through either barrier. This can be seen by assuming 
radiative boundary conditions at either infinity and calculating the tunneling probability in 
the same way as done for barrier penetration problems. As a result, there could be states 
found within the well which are quasibound. More precisely, these would be solutions with 
complex eigenvalues - the imaginary part being related to the lifetime of the quasibound 
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state in the well. Some recent methods for the calculation of the lifetime of the quasibound 
states have been proposed in We have plotted the potential for different values of the 
parameter v. It is apparent from Fig. (1) that the depth of the potential well increases with 
increasing values of the parameter ly (02). As a consequence, the probability of having a large 
number of bound states within the well gets more and more pronounced with increasing i^. 
Case II : ( tti = and a2 positive ) 

The potential obtained in this case is the well-known Posch Teller potential (Fig. 2). 
The solutions of the Schrodinger equation reveals that the energy spectrum is continuous 
for the positive eigen values whereas it is discrete for negative eigenvalues of the energy Q]. 
Discrete bound states are found inside the well. 
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FIG. 2: The potential reduces to the well known Posch- Teller potential for ai = and 02 positive. 

Case III & IV : ( ai, 02 both negative / ai -ve, 02 +ve ) 

For the single barrier potential (Case III) we essentially have the usual barrier penetration 
problem. In the last case, i.e. for a potential well, we can obtain bound states with a 
discrete energy spectrum. We intend to focus on the double-barrier volcano potential in the 
remaining part of this article. 

Before we go over to solving the Schrodinger equation let us note the following fact. The 
volcano potential in the Case I above is not entirely unrelated to the inverted double well 
which we discussed at the beginning of this section. To see this, consider a Taylor expansion 
of the terms in the cosh-sech volcano and retain terms up to order x^. One will notice the 
reappearance of the ax^ — bx^ potential for uai < 02. Polynomial approximations for the 
other cases (II, III, IV) may also be obtained using the same expansion. 
Exact solutions of the Schrodinger Equation: For the potential given in Eqn. ([21), the 
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Schrodinger equation is of the following form 



, ^ ' + {Ai cosh^'' X + A2 sech^x + A3)^(x) = (3) 

where, Ai = 2ai,A2 = 2a2 and ^3 = 2E. It is possible to find an exact solution of the 
above equation if the relations, A2 = |(| + 1) and A3 = are satisfied. The most general 
solution to Eqn. ([3]) is a combination of two linearly independent solutions, 'iIji{x) and ip2{x), 
given by 



cos \/Ai f( cosh xYdx 
^1 (^) = A ^ r I ^^ (4) 



sm 



(cosh x) 
^Ai J {cosh xydx 



M^) = B ' , (5) 

(cosh X) 2 

where, A and B are arbitrary constants, to be determined by the normalization condition. 
It is manifest from the general expressions that u should be positive for the wave function 
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to be finite at x ^ ±00. The wave function is normalizable with A = (^|(1 + e ^^^^jj 2 
and B = — e~^^^)^ ^ for u = 1 with the normalization constant chosen as unity. The 
total energy of the particle is negative and the motion of the particle is bounded between 
the classical turning points of the potential when the particle's energy lies between Vmin 
and the peak value of the potential i.e. Vmin < E < Vmax- The fact that the value of the 
bound state energy be within the well imposes a restriction of the allowed domain of Ai. 
This turns out to be > |. Let us now explore how the states depend on different values 
of u. We show this graphically. The number of nodes of the wave function increases with 
increasing u within a given region. In Fig. (3), we have plotted the wave functions for z/ = 1 
and ly = 4. The normalization constants are obtained for v = 1 and one can argue that the 
wave functions are also normalizable for v = Am. the sense that 1-01^ covers a finite area for 
X ±00. 

The probability density of the even states for the same set of v values have been depicted 
in Fig. (4). From this we can conclude that the wave functions are normalizable and well 
behaved. The nature of the wave functions and probability density for the odd states with 
different values of the parameter v are shown in Fig. (5) and (6) respectively. 

The energy of the particle states can be calculated for definite values of v from the relation, 

2 2 

A3 = — X5 which in turn implies that, = — Energy eigenvalues are negative in this 
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FIG. 3: The even parity wave function, V'i(x) is plotted as a function of x for two different values 
of the parameter = 1 and 4. It is clear that the number of nodes increases with v increasing. 
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FIG. 4: The probability density is plotted as a function of x for two different values of the parameter 
V = \ and 4 for even states. As both the curves are showing finite area we can conclude the wave 
functions are normalizable. 

case. The even and odd parity states are degenerate though the model is one dimensional. 
This is an unusual feature which contradicts the usual theorems on definite parity bound 
states in symmetric potentials. We have not been able to obtain a better understanding of 
this curious feature of our solutions. However, the singularity in the potentials may provide 
the key to the understanding of degeneracy in one dimensions [9']. 

Furthermore, one can also evaluate the expectation values of x and p (using standard 
quantum mechanics) for these states. It turns out that: 

(a;)=0 ; (p)=0 (6) 
However, one can check that (x^) is finite, whereas (p^) diverges to positive infinity. This 
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FIG. 5: The odd parity wave function, ■02 (x) is plotted as a function of x for two different values 
of the parameter v = 1 and 4. It is clear that the number of nodes increases with u increasing. 
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FIG. 6: The probability density is plotted as a function of x for two different values of the parameter 
u = 1 and 4 for odd states. As both the curves are showing finite area we can conclude the wave 
functions are normalizable. 

implies that the uncertainty Ax is finite, but the uncertainty Ap is infinite. The curious 
fact is that, the kinetic energy T has a positively infinite expectation value, while (V) is 
negatively infinite but, (H) is finite. This happens because, in the integrand, the contribu- 
tions from T and V which make their individual expectation values diverge, have opposite 
signs and therefore, cancel! 

What can we say about a bound state which is infinitely spread out in position space? 
One notices that the even though the wave function is damped the number of nodes of the 
odd and even wave functions, within a small interval dx increases with x. More specifically, 
we have 
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fAij {coshxYdx (7) 
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for the odd and even states, respectively. One can check easily that ^ becomes progressively 
larger for large x and fixed dx. This could be a reason for the divergence of the (T). In 
other words, the momentum space wavefunction has a significant contribution from the large 
momentum modes even though we have a localised bound state in position space [1C|. It 
may be noted that there are positive energy bound states (in the continuum), the so-called 
von- Neumann- Wigner (vNW) states, which do have infinite number of nodes However, 
apart from being positive energy states, note that for the vNW states, the potentials are 
bounded below unlike the potentials we have studied here. A better understanding of the 
peculiar bound states discussed in this article is surely desirable. 

Concluding remarks : Let us now summarize our results. We have provided prescriptions for 
constructing volcano potentials which seem to have, of late, gained relevance in various areas 
of physics. After a brief discussion on classical mechanics in a volcano potential we have 
studied the quantum features in some detail. For a one-parameter family of such potentials 
we have been able to write down exact bound state solutions for specific values of the energy 
and the parameters in the potential. These solutions are normalizable in the usual sense 
(square integrable) though we find degenerate states of even and odd parity. We have shown 
through plots, the nature of the wave functions and probabilities for several cases. We have 
also calculated expectation values of position and momentum, the uncertainties in x and p 
and the density of nodes for our wavef unctions. 

A point to note is that the potentials we have discussed lead to non-Hermitian Hamilto- 
nians (due to the divergence to negative infinity of the potential). Though still a matter of 
debate, a fair amount of current work has shown some of the distinctive characteristics (eg. 
real eigenvalues for complex Hamiltonians. etc.) of such Hamiltonians (the so-called PT 
symmetric theory Gj). The definition of observables as well as the construction of quantum 
mechanics needs modifications in the context of non-Hermit icity. Moreover, their physical 
utility remains an open question though we have mentioned some of them here. In fact, de- 
spite repeated use in high energy physics (braneworld models) the aspect of non-Hermiticity 
has never been mentioned in any article. 



We hope the results derived in this paper may be useful in some way to researchers 
interested in pursuing the physics of volcano potentials. 

SK thanks R. Parwani for discussions and some useful insights on various issues related 
to this article. RK thanks CSIR, India for financial support through a Senior Research 
Fellowship. 
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